Ballistic quantum spin Hall state and enhanced edge backscattering 

in strong magnetic fields 
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The quantum spin Hall (QSH) state, observed in a zero magnetic field in HgTe quantum wells, 
respects the time-reversal symmetry and is distinct from quantum Hall (QH) states. We show 
that the QSH state persists in strong quantizing fields and is identified by counter-propagating 
(helical) edge channels with nonlinear dispersion inside the band gap. If the Fermi level is shifted 
into the Landau-quantized conduction or valence band, we find a transition between the QSH and 
QH regimes. Near the transition the longitudinal conductance of the helical channels is strongly 
suppressed due to the combined effect of the spectrum nonlinearity and enhanced backscattering. It 
shows a power-law decay B~ 2N with magnetic field B, determined by the number of backscatterers 
on the edge, N. This suggests a rather simple and practical way to probe the quality of recently 
realized quasiballistic QSH devices using magnetoresistance measurements. 

PACS numbers: 



Introduction.- Recently novel two-dimensional (2D) 
electronic state - quantum spin Hall (QSH) state - has 
been theoretically proposed [l-Q and experimentally re- 
alized in HgTe quantum wells (QWs) |^-|6|- It originates 
from spin-orbit band splitting and is characterized by 
time-reversal invariant gapless states on sample edges, 
where electrons with opposite spins counter-propagate, 
while the bulk states are fully gapped. Such (helical) 
edge channels make the QSH insulators topologically dis- 
tinct from ordinary band insulators, and hold promise 
for reversible manipulation of spin-dependent quantum 
transport. Similar edge states and transport have been 
discussed in the quantum Hall regime in graphene fj- 11|. 

In experiments on HgTe QWs [3, |i| , the QSH regime 
was detected by measuring the longitudinal conductance 
of two spin channels propagating in the same direction on 
opposite edges of the sample. This finding was further 
substantiated by the observed suppression of the edge 
transport in a magnetic field 0, H| , which breaks the 
time-reversal symmetry of the QSH state, thus revealing 
the helical edge channels. The magnetoresistance mea- 
surements 0, ll[ and theory 12 1 have so far been done 
for large disordered samples. In view of the progress in 
the miniaturization of the QSH devices @ there is an ap- 
parent need to investigate the magnetotransport in the 
ballistic QSH regime, which is the goal of our work. 

The interest in the ballistic QSH transport originates 
from the profound difference between the helical edge 
states, characterized by the dissipative longitudinal con- 
ductance ~ 2e 2 /h 04a], and dissipationless chiral quan- 
tum Hall (QH) channels [HQ. As a new test to demon- 
strate this distinction, we propose to measure the longitu- 
dinal magnetoconductance of a ballistic HgTe QW when 
its edge spectrum changes from helical to chiral. Such 
a transition is expected when the Fermi level is driven 
from the band gap into the Landau-quantized conduc- 
tion or valence band where a dissipationless QH state 



sets in. The latter is insensitive to the edge backscatter- 
ing 13], Il4| . while on the QSH side of the transition we 
find strong suppression of the two-terminal conductance 
g due to the backscattering in a magnetic field B: 



g(e,B)cx(E g -\e\f N /B 



} 2N 



E n 



(1) 



Here e indicates the position of the Fermi level with re- 
spect to the middle of the band gap (equal to 2E g ). 
This result contrasts the zero-field conductance which 
increases as the Fermi energy is pushed into the metallic- 
type conduction or valence band 041] ■ Also, unlike the 
exponential B decay in strongly disordered systems [l2| , 
Eq. ([T]) describes a power-law magnetoconductance. 

Equation ([]} assumes the presence of a few (N) 
backscattering centers (see, also, Fig. [1]), such as sample 
inhomogeneities where electronic trap states can interact 
with the edge channels randomizing their propagation 
directions [6J. Although in a zero field this effect is be- 
lieved to be weak, we show that near the QSH-QH tran- 
sition the backscattering is dramatically enhanced due to 
the reduction of the group velocities of the coupled QSH 
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FIG. 1: (Color online) Two-terminal QSH device in a strong 
magnetic field. Edge channels are localized within the mag- 
netic length A. Gray regions schematically indicate backscat- 
tering centers (see, also, Eq. (|14p ). We assume that the cur- 
rent is carried by the right moving states (solid arrows ±), 
populated in contact L and equilibrating in contact R. 



2 



modes. According to Eq. ([T]), the analysis of the power 
of the magnetoconductance decay can be a simple tool 
to determine the quality of the QSH devices, which is an 
important practical task. 

Model- We will first analyze the edge states in 
scattering-free HgTe QWs using the effective 4-band 
model derived in Refs. 0,0]. In this approach one works 
in the basis of the four states near the r (k = 0) point of 
the Brillouin zone: |ei+), \hi+), |ei — ) , and [hi — ), where 
ei and hi are the s-like electron and p-like hole QW sub- 
bands, respectively. The index r = ± accounts for the 
spin degree of freedom. The effective two-dimensional 
Hamiltonian can be approximated by a diagonal matrix 
in r space 0, 0|: 



H 



h u 
h 



-k 



, hk = d k er, d k = (hvk x , -hvky, M).(2) 



where Pauli matrices o~ x ^ y ^ z act in subband space, v ~ 
5.5 x 10 5 ms~ 1 is the effective velocity and M deter- 
mines the band gap E g = \M\ at k = 0. In Eq. ([2]) we 
omit terms cx k 2 which are small near the T point and in 
the range of fields we consider [l5j]. We also neglect the 
bulk inversion asymmetry because we will focus on strong 
magnetic fields where the subband mixing is suppressed. 
Up to a unitary transformation, Eq. ^ is equivalent to 
a massive Dirac Hamiltonian Hjj — Kvt z &\s,-\- Mt z o~ z [t z 
is the Pauli matrix in spin space] . We will work with the 
corresponding retarded Green's function defined by [el— 
H d ]G(t,t') = I6(r - r'), where k = -iV - eA(r)/cft, 
A(r) = {—By, 0, 0) is the vector potential of an external 
magnetic field B, and I = r ao = diag(1, 1, 1, 1). As- 
suming a sufficiently wide sample, we find G(r,r') near 
one of the edges, e.g. y = 0, using the boundary condi- 
tion G(r,r')\ y=0 = r a x G(r, r')\ y=0 , equivalent to con- 
finement by infinite "mass" at y < [Hj], The 
matrix G — diag(G+, GL) is diagonal in r space, and 
each G T can be diagonalized in e,h space: 



G T 



re+AI 




(3) 



Expanding G in plane waves e lkx yields the boundary 
problem for the diagonal elements: 



dl- L 
d z G 



z k f 



G 



ee\rk 



A( € +tM) 

h 2 v' 2 



6{z- 



- aC , I n - A ( T£ + M ) 

ee\rk H ^-*"ee|rfc | z —q ' " ft v 



Xk, 



(4) 
(5) 



with z = y/X, Zk = — 2A k sgn(eB), A = ^ch/2\eB\, and 
a = X 2 (M 2 - e 2 )/h 2 v 2 - sgn{eB)/2. For G hh \ Tk one 
replaces r, k, B — > —t, —k, —B. The solution is obtained 
in terms of the parabolic cylinder function U(a,z) [22j |. 
It contains the edge contribution of the following form: 



a(z, z') 



G T k 



1 P(z') 



tM - T(hv/X)U(a, -z k )/U(a + 1, -z k ) ' 



(G) 



(a) M > 7 EAM\ 

x=+r 

Insulator 



T=-L 



(b)M<0 ,6/lMI 
QSH 



Zk 




Zk 





1.5 


' \£/IMI 


(d) 






- - chiral 


1 
















- helical 


0.5 






Z 






- 


Z+ 2 Z- 


3 4 



FIG. 2: (Color online) Edge-state energy vs. center-of- 
oscillator coordinate zt — — 2Afcsgn(eB) [ftw/A|M| = 0.5, 
eB > 0]: (a) insulator and (b) QSH state with two coun- 
terpropagating gapless spin channels, (c) Spatial behavior of 
zero-energy states in weak (dashed curves) and strong (solid 
curves) fields [see, Eq. ([6]) and text], (d) Numerical (points) 
and analytical (Eq. 1111 solid lines) results for the spectrum 
in strong fields, hv/X\M\ > 10. It transforms from helical to 
chiral at e = \M\. 



_ U(a,z-z k )U(a,z' -z k ) a _ U(a,-z k )U(a+l.z~z k ) 
~ \U(a,-z k )U(a+l-z k )'P U(a+l,-z k )U(a,z-z k ) ' 



(7) 



The new feature of this solution is that it is valid for an 
arbitrary parameter hv/X\M\ which measures the mag- 
netic field strength. Below we compare weak- and strong- 
Held regimes defined by hv/X\M\ < 1 and hv/X\M\ 1. 

Weak- vs. strong-field QSH channels.- The edge-state 
spectrum is given by the pole of Eq. ((6]). For weak 
fields we reproduce the transition from the band insu- 
lator with M > to the QSH state with M < [cf. 
Figs.[5Ja) and (b)], which is observed at the critical QW 
thickness sa 6.3 nm 0, IH [B|. The QSH state has two 
gapless counter-propaga ting spin modes exponentially lo- 
calized at the edge [2|, ll7(, as seen from Eq. ([6]) and 
Fig. c) where we use the asymptotic formula U(a, z) 
V^/[2 Q / 2+1 / 4 r(3/4 + a/2)]e J ^ z with |a| > 1 [18], valid 
for low fields and energies |e| < |M|: 



G T k 



(ao + ^)Me-|M|fe +3 /')/ fi , hv 



e - tMQ(M) - rhv{k - k B ) ' A|M| 



«1. (8) 



The subgap dispersion is linear: e rk — rhv{k — ks) for 
M < 0. The magnetic field only shifts the zero-energy 
point ks = — eBv/(2c\M\) with no effect on transport. 

As the magnetic field does not open a gap, the QSH 
state persists in strong fields 7w;/A|M| » 1, though 
the QSH channels are no longer localized at the edge 
[see, solid curves in Fig. [2fc)]. The electron function 
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G ee (z, z) oc a(z, z) for eB > [or the hole one Ghh(z, z) 
for eB < 0] behaves almost like the lowest-Landau- 
level bulk wave function peaked at the center of os- 
cillator (COS) 2^. The other functions are small at 
z ~ Zfe. The strong- field asymptotic is obtained for 
U(a, z) w U(-l/2, z) = e" z2 / 4 , U(o+l, z) Z7(l/2, z) = 
e z2 / 4 yV2erfc(z/V2), and /3 < 1 in Eqs. © and ©: 

Cr r /c « G T fe, G T fc(z,zj = , (9) 



a(z, z') p 
e r/c = tM + r 



2 e" 



" Aerfc(-z fe /v / 2) 
2~ fru e" 2 '/ 2 



7T A erfc(-a*/V2)' A|Af| 



(10) 



» 1, (11) 



where erfc(z) is the complementary error function. How- 
ever, the most essential distinction of this regime is the 
nonlinear spectrum (jTTJ) . Upon crossing the gap energy 
Eg = \M\ it changes from helical to chiral, as illustrated 
in Fig. H[d). Therefore, the QSH state transforms into 
a dissipationless v = 1 QH state ID, 14 1. Unlike related 
work on HgTe QWs @, [H, (H and graphene @, El, 
we intend to study the QSH-QH transition in the en- 
ergy (e.g. gate voltage) dependence of the longitudinal 
conductance. For that purpose, we need the group ve- 
locities, v±(e,B) and COS coordinates, z±(e,B), which 
are obtained from Eq. (fTT|) linearized near given energy, 
e T k ~ e — (hv T /2\)(zk — z T ). Here z T (e, B) is the solution 
of equation e T k = e, which is related to the velocity by 



v T (e, B) = 2X(t\M\ + e) z T (e, B)/h. 



(12) 



The edge state can be described by the one- 
dimensional Green's function, G T (x, x') = j e lk i x ~ x ) 
x Jg 00 dy G T k(y,y), where G T k{y,y) is localized within A 
[see, Eq. (|10p ]. Using the linearized dispersion we find 

pik T (e,B)(x — x') 

zn|w T (e,ij)| 2A 

The step function Q([x — x'\t) accounts for the chirality. 

Edge backs cattering and magnetoconductance.- We now 
calculate the two-terminal conductance of a QSH system 
[see, Fig. Q] using the scattering matrix formalism. Since 
the edges are assumed decoupled, it is sufficient to do 
the calculations for one of them, e.g., for the lower edge 
in Fig. Q] which is described by the following S matrix: 

5 = C++ 1+- ) ® ££ 4 £i . Here r's and i's are the reflec- 

RL R.Ft 



tion and transmission amplitudes for the right ("+")- 
and left ("-") -moving states; (ctq + cr z )/2 projects the S 
matrix on the electron QW subband which has the non- 
vanishing wave function [see, Eq. ([9])]. The conductance, 
g = (e 2 /h) \t^\ , is calculated using Fisher-Lee rela- 
tion 0, t++ = ih\v + \g ++ {x G R,x' e L), between <++ 



and the diagonal element Q ++ (x, x') of the Green's func- 
tion G(x, x') = ( g ++ g + ~ ). Its off-diagonal part is due to 
backscattering. We model it by the sum of N potentials, 
localized at positions x n with non-zero matrix elements 
V n between the right- and left-moving states: 



L — J n—1..N 



(14) 



Microscopically, the coupling between the counter- 
propagating channels can be mediated by interaction 
with electronic trap states which are likely to exist even 
in high quality samples [6]. Note that choosing the 
other off-diagonal matrix, r y does not change the final 
result. Potential (fT4f results in the Dyson equation 

G(X,X') = G(X,X') + Y.n=l..N G ( X ^ X n)V n T x G(x n ,x'). 

This allows us to express the off-diagonal part 
Q T .- T (x, x') through Q TT {x,x r ) and obtain a closed 
equation for the latter: Q TT (^x,x / ) = G T (x,x / ) + 

Z^n m—l..N ^ T ' Xn ) VnG— T (x n , X m ) VrnQrr {Xm ) X ) . 

With known unperturbed function G T (fT3]) and for not 
large N, we solve this equation and calculate g. Let us 
look first at the particular cases N — 1, 2 and 3: 



Yt - 
h 2 \ vjf-v— | - 



_ e 2 li , V 1 2 + Vi+ViV2c' Qd12 
~ h I 1 i K 2 \v+v-\ 



(15) 
(16) 



£ V 2 +V!V 2 C iQd ^ + V 1 y 3 c iQd 13 +V 2 V 3 O iQd 23 



~> w\v + v-\ 
y 2 v£ +y 2 y 3 2 + y 2 2 y 3 2 +y 2 y 2 v 3 c iQd23 , y 2 y 2 2 y 3 2 

H ft 4 U^'t>^ fi.6|„ +aJ _ |3 



,(17) 



where Q = k + — k- and d nm = x m — x n . 

It is clear from Eqs. (|15[) - pT|) that for an arbi- 
trary TV the conductance contains the cross product 
Vi ■ ■ ■ Vh/\v+V-\ n arising from the simultaneous scat- 
tering from N potentials. This is the most divergent 
term when one of the velocities v± vanishes near the band 
gap, |e| -> E g = \M\ (e.g. u_ -> in Fig. Hfd)). Such 
strong enhancement of the backscattering leads to the 
suppressed conductance, 



g » (e7/i) x ft 4iv | u+u _|^ v /(Ui • • • Uat) 4 « e'/h. (18) 

Using Eq. (|12p for and ommiting the slower func- 
tions z±(e, B), we obtain the qualitative energy and field 
dependence of the conductance near the QSH-QH tran- 
sition, presented in the introduction [see, Eq. ((TJ]. 

In a wider range of energies and fields the typical 
behavior of the conductance can be understood from 
Eq. (fTB]) assuming two backscattering centers on the 
edge. First of all, it is easy to verify that Eq. (fl6|) is 
valid not only for strong fields, but also in the weak-field 
case where the unperturbed Green's function is given by 
Eq. (JSJ). Since the weak- field spectrum is linear e T k — 
rhv(k + ks), we have v + = —v- = v, k± = —ks ± e/hv 
and Q = k + — fe_ = 2e/hv. Therefore, g is independent 
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FIG. 3: (Color online) Conductance [see, Eq. (|16|l ] vs. (a) 
energy below the band gap E g = 1.5 meV and (b,c) vs. mag- 
netic field. Vi,2 are in units of meV-^m; d\2 = 3/xm. 

of the magnetic field and for Vj. 2 "C is almost inde- 
pendent of energy (see dashed curve in Fig. 3a). Thus 
weak channel mixing is hardly detectable for small B. 
In contrast, in strong magnetic fields, scattering of the 
same strength is sufficient to suppress the conductance, g 
near the band gap where the QSH-QH transition occurs 
[cf. dashed and solid curves for V1.2 = 0.01 meV-/zm 
near E g = 1.5 meV in Fig. [3ja)]. Figs [5Ja),(b) and 
(c) also show that the conductance suppression is ac- 
companied by Fabry-Perot-type oscillations due to in- 
terference of the counter-propagating channels which ac- 
quire the energy- and field-dependent phase difference 
(fc + — fc_)rfi2, in scattering between the defects. 

In our model the upper magnetic field limit lies in 
the range of a few Tesla. This estimate is based on 
Refs. 0, [TU predicting another (B-field induced) QSH- 
QH transition due to a quadratic correction 2?k 2 to the 
mass term in effective Hamiltonian J2]). The smallness 
of the parameter B\M\/2h 2 v 2 <C 1 [a El allows us to 
neglect such k 2 term and to meet, at the same time, the 
strong field condition hv/X\M\ > 1. 

Conclusions.- We have studied the longitudinal con- 
ductance of helical spin edge channels in HgTe quantum 
wells in strong magnetic fields. As the Fermi level ap- 
proaches the band gap, the conductance vanishes due 



to the combined effect of the spectral nonlinearity and 
channel backscattering. Such energy dependence indi- 
cates a transition between the quantum spin Hall and 
dissipationless quantum Hall regimes. The conductance 
exhibits B~ 2N magnetic field dependence, determined by 
the number N of backscattering centers on the edge. This 
suggests a simple way to detect defects in ballistic QSH 
devices using standard magnetoresistance measurements. 
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Tool r , - r°° (? ?'\ r a ^ k u i a > z k)+<iu^^k) u( 
u ee \ Tk - Or ee | Tfc (z, z ) - o d ^m a> ^ Zk)+qU(a ^ Zk) u z 

Zk)U(a,z' — Zk). The last term is the edge state, while 

G^, Tfe = C[Q(z - z')U(a, z - z k )U(a, -z' + z k ) + Q{z' - 

z)U(a, z' — z k )U(a, — z+Zk)] is the bulk solution with C = 

-A(e + rM)P(a+ l/2)/ v / 2^fi 2 u 2 . For eB > and using 

recurrence relations for U(a,z) we obtain Eq. ((6j). 



